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We analyze the separability properties of density operators 
supported on <E 2 ® (P whose partial transposes are positive 
operators. We show that if the rank of p equals N then it is 
separable, and that bound entangled states have rank larger 
than N. We also give a separability criterion for a generic 
density operator such that the sum of its rank and the one of 
its partial transpose does not exceed 3N. If it exceeds this 
number we show that one can subtract product vectors until 
decreasing it to 3N, while keeping the positivity of p and 
its partial transpose. This automatically gives us a sufficient 
criterion for separability for general density operators. We 
also prove that all density operators that remain invariant 
after partial transposition with respect to the first system are 
separable. 

I. INTRODUCTION 

Perhaps, entanglement is the most intriguing property 
of Quantum Mechanics. It arises when two or more sys- 
tems are in a non-separable state; that is, in a state that 
cannot be prepared locally jlj]. Apart from having fun- 
damental implications ||, such as Bell's theorem or the 
existence of decoherence, entanglement is in the realm 
of several practical applications of Quantum Information 

Despite of its importance, we do not know yet how to 
quantify entanglement. In fact, even in the simple case 
in which we only have two systems A and B in a mixed 
state p, there exists no general criterion that allows us 
to distinguish whether the state is separable or not; that 
is, whether it can be written as a mixture of product 
vectors of the form |e, /) = \e)A <8> \ J)b- Several nec- 
essary conditions for separability have been formulated 
in the recent years. Werner jjj formulated a criterion 
based on the analysis of a local hidden variable model 
and on the positivity of the mean value of the flipping 
operator. The Horodeckis obtained a criterium in the 
form of inequalities for the so called a-entropies. An im- 
portant step forward to clarify the situation was taken 
by Peres ||. Peres criterion is based on the positivity of 
the operator obtained after partially transposing the den- 
sity operator. This criterion gives a necessary condition 
for separability for density operators in Hilbert spaces of 
arbitrary dimension: if the partial transpose of a given 
density operator p is not positive, then p is not separable. 

The Horodecki family || has established the connec- 
tion between separability and the theory of positive op- 



erator maps Jj . A necessary and sufficient condition for 
separability has been expressed in terms of these maps; 
such condition, unfortunately, is constructive only if the 
structure of the set of positive maps is sufficiently known. 
This turns out to be the case for the situations in which 
the Hilbert space corresponding to the first subsystem is 
two-dimensional and the one of the second is either two- 
or three- dimensional, the so-called 2x2 (two qubits) and 
2x3 cases (one qubit and one qutrit), respectively. Using 
that result, the Horodeckis have shown that Peres crite- 
rion (called hereafter Peres-Horodecki criterion) provides 
a sufficient condition for separability in 2 x 2 and 2x3 
systems. The fact that this is not a sufficient condition 
for higher dimensional systems was pointed out later on 
||, by giving examples of density matrices in 2 x 4 and 
3x3 spaces that are entangled, but have positive par- 
tial transposes (PPT). After that, a systematic method of 
constructing entangled states with PPT, based on the ex- 
istence of unextendible basis of product vectors was for- 
mulated [3,[L0|. Recently, positive map theory has been 
used to obtain a reduction criterion of separability pd| ] . 

So far, there have been very few results concerning 
sufficient conditions for separability in high dimensional 
systems |T^,|l3|]. In this paper we analyze the separabil- 
ity properties of density operators acting on <D 2 ®(E N . 
This situation is particularly useful since it allows to de- 
termine a number of separability properties in the mul- 
tiqubit case |L4|]. Our analysis is based on subtracting 
product vectors from p while keeping p and its partial 
transpose positive. The idea of subtracting product vec- 
tors was developed in |]l5|] and has proven to be very 
useful so far since it allowed: i) to define unique optimal 
decomposition of mixed states in Hilbert space of arbi- 
trary dimensions into a sum of separable and insepara- 
ble part; ii) to find minimal decompositions into product 
vectors of separable density operators acting on (D 2 <S> <C 2 
(two qubits), and iii) to define a measure of entangle- 
ment in those systems. In the present paper, by studying 
the conditions under which we can subtract such product 
vectors we are able to find necessary and sufficient con- 
ditions for separability for low and medium rank density 
operators. We also give a sufficient condition for sepa- 
rability for general density operators. Furthermore, we 
show that when quantum systems remain invariant after 
partial transposition with respect to the first system then 
they are separable. This last result displays the differ- 
ence between the 2 x N case and the general M x TV one 
(N,M > 2), where it has been found that there are op- 
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erators satisfying this condition which are not separable 
§• 

The structure of this paper is the following: In Sec- 
tion II we give the definitions that we use throughout 
the paper. In Section III we derive a series of lemmas 
concerning product vectors. First, we recover the results 
given m Ref. II concerning the conditions under which 
one can subtract product vectors from a density opera- 
tor. We then give a separability criterion for the case in 
which there is a finite set of product vectors that can be 
subtracted. Next we show that if there are product vec- 
tors in the kernel of p, then the problem of separability 
can be highly simplified. Finally, we analyze under which 
conditions there exist product vectors in subspaces (e.g. 
the kernel and the range of p), and show that the prob- 
lem reduces to finding the roots of complex polynomials. 
Section IV contains most of the results of this work con- 
cerning separability. We first show that one just have 
to consider the cases in which the rank of p is larger or 
equal than N. When the rank of p is equal to N (and 
p is not trivially supported on a smaller subsystem of 
B) then it is separable. From this result follows that 
bound entangled states have a rank larger than TV. We 
continue by analyzing the case in which the sum of the 
ranks of p and its partial transpose is smaller than 37V, 
giving a necessary and sufficient separability criterion for 
generic density operators and for N < 10. After that, we 
show how in the general case one can reduce the problem 
to that situation. Then, we prove that if p is equal to 
its partial transpose with respect to the first subsystem, 
then it is separable. The same statement holds when p is 
"sufficiently close" to its partial transpose (in the sense 
of operator norm, for instance) . Finally, in Section V we 
illustrate our results by analyzing the case (D 2 <8> <C ■ 

II. SEPARABILITY IN 2 x N 

We consider a density operator p > acting on (D 2 <S> 
(U N . We wish to find out conditions such that we can 
state whether p is separable or not; that is, whether it 
can be written as 

P = ^2\ei,fi){ei,fi\, (1) 

i 

where \ei, fi) = \ei) A ® \fi)B are unnormalized product 
states. Throughout this paper we will make use of the 
fact that a necessary condition for separability is that 
the partial transposition of p is a positive operator [^J. 
Given an orthonormal basis {|0}, |1}} S (f 2 , we define the 
partial transposition of p with respect to the first system 
in that basis as 

l 

P Ta = E Mi\p\j)A\j) A (i\ (2) 

i,3=0 



Thus, from now on we will assume that both p, p A > 
(otherwise p is already not separable). 

Throughout this paper K(X),R(X),k(X), and r(X) 
denote the kernel, range, dimension of the kernel, and 
rank of the operator X, respectively. By |e*} we will 
denote the complex conjugated vector of |e) in the basis 
0}^, \1}a in which we perform the partial transposition; 
that is, if |e) = a\0)+/3\l) then |e*) = a*\0) + /3*\l). The 
subscript "r" will be used to denote real vectors; that is 
|e r ) = |e*). We will denote by |e) the vector which is 
orthogonal to |e). On the other hand, we will make use 
of the following property, 

A(e i \p\e j ) A = A (e*\p T -\e*) A (3) 

for any pair | e^) , \ej) £ (T 2 , which can be easily deduced 
from the definition of partial transposition. Note that if 
p is separable [cf (Q)], then 

/^ = Ei e ^>«</*i' ( 4 ) 

i 

is also separable; obviously, the opposite is also true. 

We say that p acting on <E 2 ® (E is supported in 
(T 2 ® <U M if the minimal subspace H C (D N such that 
R(p) C (C 2 ® H has dimension M. This means that there 
exist N — M (and not more) linearly independent vectors 
{\fi),i=l,2,...,N-M}e(D N such that \e,fi) € K{p) 
for all |e) £ (E 2 . The utility of this definition can be seen 
in the following example. Let us consider an operator p 
acting on (D 2 (g)<Z7 2 . We can always trivially regard this 
operator as acting on a higher dimensional space (e.g., 
(D 2 ® <F 4 ) but giving zero whenever it acts outside the 
original space. If we are studying the separability prop- 
erties of p it is useful to get rid of the spurious dimensions 
and really consider p as acting on (D 2 ®(D 2 . Thus, we will 
concentrate on operators not only acting, but also sup- 
ported on <E 2 ®<E N . 

The method we use to find separability conditions is 
based on subtracting product vectors from p in such a 
way that both p and p TA remain positive. If we are 
able to subtract several of those product vectors until 
we achieve that the result vanishes, then we will have 
obtained that p is separable. Thus, a crucial point is to 
find the conditions under which we can find such product 
vectors. In the next few sections we will determine those 
conditions. 

III. PRODUCT VECTORS 

As mentioned above, our procedure to study separa- 
bility is based on subtracting product vectors |e, /) from 
the density operator p in such a way that the resulting 
operator and its partial transpose remain positive. Thus, 
we have to determine under which conditions this is pos- 
sible. This section deals with that problem. We have 



2 



divided it in three parts: (i) Properties related to the 
existence of product vectors in the range of p or p TA ; 
(ii) Properties related to the existence of product vec- 
tors in the kernel of p or p T - A ; (iii) Existence of product 
vectors in subspaces. In the first subsection we derive 
some useful properties of a density operator whenever 
there are product vectors in its range. In particular, we 
will recover the results of Ref. |L5| which state that we 
can only subtract product vectors of p that belong to its 
range, and that by subtracting the right amount we can 
always decrease its rank. Then we will show that if we 
have a set of product vectors in the range of p and an- 
other related set in the range of p TA , then p is separable 
iff it can be written as a mixture of those vectors. This 
means that if we can find all the product vectors satis- 
fying these conditions, then we can check whether p is 
separable or not by using only those vectors. In the sec- 
ond subsection we will derive some other properties of p 
whenever there are product vectors in its kernel. In par- 
ticular we will show that in that case one can subtract 
one product vector and not only reduce the rank of p, 
but also reduce the problem to finding whether a related 
density operator supported on (D 2 x (Z^ -1 is separable. 
This means that the product vectors in the kernel play a 
very important role since they allow to simplify the prob- 
lem of separability. Finally, in the last subsection we will 
derive the conditions needed to have product vectors in 
subspaces; in particular in the kernel and the range of p 
and p TA . We will see that in order to find them, one has 
to solve some complex polynomial equations. The solu- 
tions to those equations determine which are the product 
vectors that can be subtracted from p, which, together 
with the results of the first subsection, allow us to an- 
alyze the separability of certain density operators, as it 
will be shown in the next section. 



A. Range of p 

We have divided this subsection in two: In the first one 
we introduce a lemma which gives the conditions under 
which the difference of two positive operators is a positive 
operator itself. This allows us to recover in a simple way 
the results given in fl5f| regarding subtracting product 
vectors from p. In the second one we show that the set 
of (common) product vectors in the ranges of p and p TA 
suffices to study the separability properties of p. 



1. Subtracting product vectors 

We use the operator norm, namely \\A\\ = 
maxii|0\|| = i||A|0)||. For a given hermitian operator X we 
will denote by X^ 1 its pseudoinverse; that is, if the spec- 
tral decomposition of X — Xk\xk){xk\ with Xk 7^ 0, 
thenX-^E^VfcX^I- 



Lemma 1: Given two hermitian operators X, Y > 
then X - Y > iff K{X) C K(Y) and \\Y^X~^ || 2 < 1. 

Proof: (if) We can write an arbitrary vector \^) = 
|* r ) + |* fc ), where |* r ) G R(X) and G K(X) C 

K(Y). Therefore (V\X - Y\V) = (W\X - Y\W r ). Since 
\V) G R{X) = R(X 1 / 2 ) we have that there exists some 
|$) G RiX 1 ' 2 ) such that |* r ) = X- 1 / 2 ^). Using that 



(<I>\X- 1/2 YX- 1/2 \<P) < ($|$), 



(5) 



we obtain that (^ r \X - Y\$ r ) > 0, and thus X - Y > 
0. (only if) First, K(X) C K(Y) since otherwise there 
would be a vector |*) G K{X) for which (*|Y|*) > 0, 
and therefore (^\X-Y\^) < 0. Moreover, V|* r ) G R(X) 
we can find |$) G RiX 1 / 2 ) such that |* r ) = X 1 / 2 ^). 
Using that > ($|Y~|$) we obtain that 



1 > (* r |^ _1/2 ^ _1/2 |* r )/(* r |* r ) 1 



(6) 



which immediately gives 

We can use this lemma to recover the results of Ref. 
|l5| . Given an operator p, a product vector \e, f) and 
A > 0, we define 



p = p- A|e,/)(e,/|, 

A ° = (e,/|p-i| e ,/)' 
1 



A 



(7a) 
(7b) 

(7c) 



We then have 

Corollary 1: p, p TA > iff \e,f) G R(p), \e*J) G 
R{p TA ) and A < min(Ao, Ao). 

Proof: We can simply use Lemma 1 with X = p 
and Y = A|e, /)(e, f\, and with X = p TA and Y = 
A|e*,/)(e*,/|.D 

The most interesting situation occurs when we choose 
A = min(Ao,Ao), since in that case we can reduce the 
ranks of p and/or p TA , as states the following 

Lemma 2: Let \e, f) G R{p), \e*,f) G R(p TA ). Then 

(i) If A = A < A then r(p) = r{p) - 1 and r(p TA ) = 

r(p TA ). 

(ii) If A = A < A then r(p) = r(p) and r(p TA ) = 

r(p TA ) - 1. 

(iii) If A = A = A then r(p) = r(p) - 1 and r(p TA ) = 
r(p TA ) - 1. 

Proof: Here, we just prove (i) , since for the other two 
cases one can follow the same arguments. First, note 
that r(p) > r(p) — 1 since they differ only by one one- 
dimensional projector. On the other hand, K(p) C K(p) 
for if I*) g K(p) then p\^) = -X {e,f\^>) = since 
e, /} G R(p) -L K(p). Besides, the vector /9 _1 |e, /) is not 
in K(p) but it is indeed in K{p), which can be checked 
by substitution in Eq. (f7a|). □ 
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This lemma allows us to decrease either the rank of p, 
of its partial transpose, or of both of them by subtracting 
a product vector |e, /} which belongs to R(p) and such 
that \e*,f) is in R(p TA ). 

2. Separability and product vectors 

Let us denote by {\ei, fi),i = 1,2, ... ,L} the product 
vectors such that |ej,/j) € R(p) and |e*,/j) £ R(p TA ). 
Note that we assume here that the total number of such 
vectors L is finite. We also denote by Pi = |e i5 fi)(ei,fi\. 
We have 

Lemma 3: p is separable iff it can be written as a 
convex combination of Pi. 

Proof: (if) Use the definition of separable state; (only 
if) We just have to notice that if a positive operator p can 
be written as a convex combination of certain projectors 
on the vectors \^i), then these vectors must belong to 
its range. Otherwise, there would exist a certain vector 

I**,) = m )+\K)> where i*u e r o>) and + \K) e 

K(p), so that (** |p|*i ) > 0, which contradicts the fact 
that \^>i ) £ K(p). Thus, if p is separable then it can be 
written as a convex combination of projectors on product 
vectors \&i, fi) £ R(p). Also p TA can be written as the 
same convex combination, but with \e*,fi) £ R(p TA ).0 

Let the projectors Pi be linearly independent. These 
projectors belong to the space C[R(p)] of linear operators 
acting on R(p), which has dimension r(p) 2 . On the other 
hand, the projectors P[ A will be linearly independent as 
well, since partial transposition is a linear invertible map. 
These projectors belong to the space C[R{p TA )\, which 
has dimension r{p TA ) 2 . Therefore, if the Pi are linearly 
independent we must have L < min[r(p) 2 , r(p TA ) 2 }. We 
can always complete the set Pi with other operators Pi 
(i = L+l, . . . , r{p) 2 ) to form a basis in C[R{p)\. That is a 
Hilbert space with the scalar product (A, B) = tr^B). 
Thus, we can construct the biorthogonal basis to {Pi} 
in C[R(p)]; that is, we can find a set of operators Qi £ 
L[R{p)] such that tr(QtPj) = Then we have the 

following necessary and sufficient separability condition 

Lemma 4: If Pi are linearly independent, then p is 
separable iff tr(Q\p) >0Vi<L and tr(Q\p) = Mi > L. 

Proof: Since the Pi form a basis, p can be expanded in 
a unique way as 

r(p) 2 

with Ci — tr(Qip). (if) since c% > we have that p can be 
written as a convex combination of projectors on \e%, fi), 
i.e. it is separable; (only if) If p is separable, according 
to Lemma 3 it can be written as (^) with Cj = if i > L 
and a > if i < L.D 



This lemma is important since in general, if L < 
min[r(p) 2 , r(p TA )], the projectors will be linearly inde- 
pendent. The reason is that given a set of L < D 2 
random product vectors which span a Hilbert space of 
dimension D, then the corresponding projectors are lin- 
early independent. 

B. Kernel of p 

We will show that if there is a product vector in the 
kernel of p, then one can find a positive operator acting 
on (D 2 (gjff^ -1 which has the same separability properties 
as p. This means that one can reduce the dimensional- 
ity of the Hilbert space of the second subsystem, which 
simplifies the problem. We will use later on these results 
to prove facts regarding low rank density operators. We 
start proving a simple 

Lemma 5: |e*, /} £ K( P Ta ) iff |e, /} £ K{p). 

Proof: (if) We have = {ej\p\ej). Using |), 
we have = (e*, f\p TA \e* , f) and since p TA > then 
P TA \ e * i f) — 0; (only if) One can prove it in the same 
way.D 

Lemma 6: If |e, /} G K(p) then either: 

(i) \e,f)£K(p),ox 

(ii) there exists a non-zero \g) £ (D N such that p\e, f) = 

\e,g) andp^e*,/) = \e*,g). 

Proof: According to Lemma 5, if |e, /} £ K(p) then 
\e*,f) £ K(p TA ). Therefore 

0=(e* : h\p TA \e*,f) = (e,h\p\eJ), (9) 

V\h) £ (D N . Thus we know that either p\e,f) = or 
p\e, f) — \e,gi) for some \gi) £ <D . Analogously, since 
e, /} £ K(p) we have that either p TA \e*,f) — or 
P TA \e*,f) = |e*,<72) for some \g 2 ) £ <D N . In the first 
case, using Lemma 5 we have (i). For (ii), it remains to 
be proven that \g\) = \g2)- Using the orthogonality of |e) 
and |e) we have 

ll|e)|| 2 |5i) = <e|p|e,/) = <e>^|r,/) = |||e}|| 2 |.g 2 }.n 

(10) 

The first part of this lemma simply states that p is 
supported on a smaller space (D 2 <S> <C M with M < N, 
since \e, f) £ K(p) for all |e) £ (D 2 . In other words, that 
we have spurious dimensions in the second subsystem. 
In the following, we will not consider that trivial case, 
and therefore we will assume that p is supported on the 
whole space (D 2 x (D N . In that case, Lemma 6 together 
with Lemma 2 tell us that there is a vector that can be 
subtracted from p. In fact, this gives rise to one of the 
most important results of this section: 

Lemma 7: If p is supported in (D 2 <S> (D and there 
exists a product vector in K{p) then p = p + p Sl where 
p s is a projector on a product vector and 
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(i) r{p) = r(p) - 1 and r(p TA ) = r(p TA ) - 1. 

(ii) p is supported on (F 2 ® (F^ -1 . 

(iii) p is separable iff p is separable. 

Proof: We have that |e, /) S K(p). Using the previous 
Lemma (ii) (since p is supported in (F 2 ®(D N ) we have 
that p\e,f) = \e,g) and p TA \e*,f) = |e*, 5). According 
to Corollary 1 we can subtract the product vector \e,g), 
obtaining p = p — p s where p s = A|e, g){&, g\ and 



A = A 



1 



1 



(e,g\p-i\e,g) \\\eW(g\f) 
= A . 



\g\{p TA )-^\e*,g) 



(11) 



According to Lemma 2 (iii) we have decreased the ranks 
of both p and p TA by one, which proves (i). We also have 
that |e, /), e, /) G K{p) and therefore p is supported in 
<F 2 <g>(F M with M < N. Now we show that M = N — l. If 
it was smaller, then there would exist \h) orthogonal to 
I/) such that |e, h), |e, h) G K(p). In that case we would 
have p|e, h) = and p\e, h) = c|e, 3) where c is a constant. 
If we define /') = c\f) — \h) ^ we would have that 
|e, /'), e, /') G K(p), contrary to our assumption that p 
is supported on <F 2 ®(Z7 , which proves (ii). It remains to 



be shown (iii) that p is separable iff p is too: 
(only if) If p is separable, and since p\e, f) 
always write 



P 



\ei,fi){ei,fi\ + |e)(e| 7?, 



(if) trivial; 
= we can 



(12) 



where (f\f l ) 



= and r\ is a positive operator act- 
ing on (F JV . If we impose \e,g) — p\e, f) we obtain 
1.9) = ll|6)|| 2? ?|/), and therefore \g) E R(f]). We can write 



= E \ e iJi)( e i>M + Ie)(e| (»(»/- Xo\g){g\), (13) 



P = 



so that if we show that the operator (77 — A |<7)(g|) > 
then we have that p is separable. Using ([ll]) we have 
that such operator is 



1 



V 



ll|e)|| 2 (.9l/) 



\g)(g\ = v 



1 



(g\ri- l \g) 



\g)(gl (14) 



which is positive according to Lemma l.D 

This theorem simply says that when there is a product 
vector in the kernel or p, then we can reduce the problem 
without changing the separability properties of p. Later 
on we will consider states which are invariant under par- 
tial transposition. In that case, we can particularize this 
result and obtain: 

Lemma 8: If in addition to the premises of Lemma 
7 p = p TA then we can always obtain the results of that 
lemma with p — p TA . 

Proof: We know from Lemma 5 that if there is a prod- 
uct vector \e,f) £ K(p) then \e*,f) G K( P Ta ) = K{p). 



With these two vectors we can always construct e r , /) G 
.ftf(p) where \e r ) is real (in the basis in which we take the 
partial transposition) . If we proceed with this vector as 
in the previous lemma then we obtain the desired result. 
□ 

Finally, we will show a property of the vectors con- 
tained in K(p) for the case when there is no product vec- 
tor in K{p). We will use this property later on to simplify 
the search of product vectors in R{p). Let us denote by 
{1*1), z = l...k{p)} and {\^),i = l...k(p TA )} a ba- 
sis in K(p) and K(p TA ), respectively. Then we have the 
following: 

Lemma 9: If there exists no product vector in K(p) 
then both {(e|\T/i)} and {(e*!* 2 )} are sets of linearly in- 
dependent vectors in (D N for any |e) G (F 2 . 

Proof: We start proving that for any |e) G (F 2 , the 
operator {e\p\e) is invertible. Otherwise there would exist 
\h) G (D N such that = (h\{e\p\e)\h) and since p > then 
p|e, h) = 0, which contradicts our assumptions. On the 
other hand, since {|e), |e)} is a basis in (F 2 we can write 
the vectors in the kernel of p as 

|* J 1 ) = |e,^) + |e,*f). 
Using that they are in the kernel of p we find that 
1 



(e|p|e> 



(e|p|e)|^ e ), 



(15) 



(16) 



where we made use of the fact that the operator (e|p|e) is 
invertible. Substituting this expression in (Eq), we have 



1 



(e|p|e> 



(e|p)|e,*f). 



(17) 



Since the l^f^'s are linearly independent, so must be the 
l^l^'s. The same argumentation holds for the vectors in 
the kernel of p TA . □ 



C. In Subspaces 

As we have shown above, it is important to know when 
one can be sure that there exists a product vector in ci- 
ther the range or the kernel of a density operator. In 
this section we will prove that for sufficiently large sub- 
spaces there always exist product vectors. We will denote 
by {|0}, |1}}, and {|1), . . . , \N)} orthonormal basis in (F 2 
and (D N , respectively. 

Lemma 10: Any subspace H G (F 2 Cg) (D N with 
dim(_ff) = M > N contains an infinite number of prod- 
uct vectors. If M = N it contains at least one product 
vector. 

Proof: We denote by i = 1, . . . 2N - M} a basis 

in the orthogonal complement of H . We have 



N 



m =j2[A iik \o,k)+B itk \i,k)), 



(18) 



k=l 



5 



where A and B are (2N — M) x N- matrices. We write 
the product vector \e, /} as 



|e,/) = (a|0}^ + |l) A )®^/ fc |fc) B . 



(19) 



Imposing that e, /) is orthogonal to the l^) for all i gives 
(aA* + B*)f = 0. This can be regarded as an equation 
for /. When M > N there exist non trivial solutions for 
each a, i.e. for each |e). For M = N we have non trivial 
results provided det(a^4* +B*) = 0; the determinant is a 
polynomial of iV-th degree in a, which always has roots. 



□ 



Any subspace H C (C 2 ® (E N with 



Corollary 2 

dim(_ff) = M > N contains an infinite number of prod- 
uct vectors of the form \e r , f), where \e r ) = \e*}. 

As we have seen in previous subsections, we can sub- 
tract vectors from p while keeping the positivity of both 
p and p TA provided we can find a product vector |e, /) G 
R{p) and such that \e*J) G R(p TA ). Wc can use the 
idea given in the proof of the previous lemma to find out 
under which conditions this is possible. Given any two 
subspaces Hi,H 2 G (D 2 ®(D N with dim(f/i j2 ) = Mi, 2 we 
denote by {\^]f 2 ),h,2 = 1, . . . 2JV - Mi, 2 } a basis in the 
orthogonal complement of -ffi j2 , respectively. Then we 
have 

Lemma 11: (i) If Mi + M 2 > 3iV then there exists 
an infinite number of product vectors |e,/) G Hi such 
that \e*J) G H 2 ; (ii) If Mi + M 2 < 3 A then there 
exists a product vector |e, /) G Hi such that e*, /) G H 2 
iff we can find an a such that there are at most N — 1 
linearly independent vectors among the following vectors: 



W^JO> + <^Jl>,a*<* 2 JO) 
Proof: We have 



iV 



(20) 



where A 1 ' 2 and B 1 ' 2 are (2N-Mi >2 }xN- matrices. Writ- 
ing a product vector e, /) as in (|19| ) and imposing that 
it is orthogonal to and |e*,/) is orthogonal to |\1/ 2 ) 
for all i, we arrive at the following equations: 



[a(A 1 r + (B 1 )*]f = 0, 
[a*(A 2 y + (B 2 y]f = 0, 



(21a) 
(21b) 



which can be regarded as 4 A — Ml — M 2 equations for 
/. For Mi + M 2 > 37V, there are less equations than 
variables, and thus there always exists a solution for each 
a, i.e. for each e), which proves (i). For Mi + M 2 < 37V 
there exists a nontrivial solution iff the matrix M(a, a*) 
composed of a(A 1 )* + (B 1 )* and a(A 2 )* + (B 2 )* has rank 
smaller than N, which is equivalent to the statement of 
the lemma. □ 

Let us now analyze the case in which Mi + M 2 = 3iV. 
Here, we have that the conditions of Lemma 11 are ful- 



filled iff det[M(a, a*)] — for some a. This determi- 
nant will be a polynomial P 2 n-Mi,2N-m 2 (a, a*) of de- 
gree 2N — Mi in a and of degree 2N — M 2 in a* . Thus, 
the existence of product vectors e,/) G Hi such that 
e*,/) G H 2 is directly related to the existence of roots 
of that polynomial. For Mi + M 2 < 3iV we have to 
impose that several determinants containing N rows of 
M(a,a*) vanish; again, those conditions are equivalent 
to imposing that several polynomials in a and a* have 
common roots. 

Note that for complex polynomials there is no theorem 
which allows to know how many (if any) roots exist. For 
example, the equation a*a + 1 has no roots, whereas a 2 — 
(a*) 2 = has infinitely many. However, in a generic case, 
given a polynomial P(a,a*) one can find a polynomial 
Q(a) of higher degree in a (say K) such that the roots 
of P are also roots of Q. Since we know that Q has at 
most K solutions, it means that P can have at most K 
solutions as well. For example, if we have P(a,a*) = 
a 2 — a* — 0, we can complex conjugate this equation 



and obtain P(a, a* 



(a*) — a — 0, which, upon 



substituting a* — a we obtain Q(a) ~ a — a — 0. This 
equation has four solutions (0, 1, e l27r / 3 , e* 27r / 3 ) which are 
exactly the roots of P. Note that we have been able 
to reduce the original equation to one containing a only 
given the fact that the P* ^ P. It is clear that the 
polynomials fulfilling this condition are dense in the ring 
of all polynomials, and this is the reason why we say that 
in the generic case one can find Q{a). 

In summary, in this Section we have shown that if there 
is a product vector in K(p) then we can find another op- 
erator p s supported on (D 2 ® (C^ -1 with the same sep- 
arability properties as p. In order to check whether we 
can have such product vectors, we have shown that if the 
dimension of a subspace [e.g., K(p)\ is greater or equal 
to N, then we can always find product vectors in it. On 
the other hand, if there is a product vector |e, /) G R(p) 
with |e*,/> G R(p TA ), then one can decrease the rank 
of p, p TA or both. This is always possible if the ranks 
of p and p TA are sufficiently large. Otherwise, one has 
to find roots of polynomials. We have also shown that if 
there is a finite, sufficiently small number of product vec- 
tors fulfilling this condition, then we automatically have 
a necessary and sufficient condition for separability. 



IV. RESULTS 

This section contains the main results of the paper 
concerning separability of a density operator supported 
on (D 2 (g)(F . We have divided them in four subsections. 
The first three deal with different situations depending 
on the rank of p. Actually, we just have to consider the 
cases in which r(p) > N as the following lemma states: 



Lemma 12: If p is supported on (D z <g>(D N then r(p) > 



N. 
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Proof: Let us assume that r(p) < N. Then K(p) has 
dimension larger than N , and according to Lemma 10 it 
contains a product vector. Using Lemma 7 we can find 
an operator pi supported on(F 2 (g)(Z7 Ar_1 and with r(pi) < 
N — 1. Again, K (pi) has dimension larger than N — 1, 
contains a product vector, and we can find an operator 
pi supported on <E 2 ®(D N ~ 2 and with r(p 2 ) < N - 2. Wc 
can proceed in the same way until we have an operator 
Pn-i supported on (P 2 gifF 1 and with r(p 2 ) < 1, which is 
impossible. □ 

The first subsection shows that if r(p) = N then p 
is separable. From that follows that bound entangled 
states have rank larger than N. The second subsection 
deals with operators p such that r{p) + r(p TA ) < 3N. 
It is shown that in the generic case and for N < 10 
one can perform a constructive separability check. The 
third subsection shows that given a density operator, one 
can always subtract vectors until the condition r(p) + 
r(p TA ) < 3iV is obtained. In the fourth one, we show 
that if p — p TA then p is separable. We also show that if 
p and p TA are "sufficiently close" then p is also separable. 

A. rank(p) = N 

Theorem 1: If p is supported on (D 2 ®(D N and r(p) = 
N then p is separable. 

Proof: We use induction. For N = 1 the statement 
is true. Let us assume that it is true for N — 1. Wc 
take now p supported on <D 2 ®<F JV , with r(p) = N. Then 
K(p) has dimension N, and according to Lemma 10 it 
contains a product vector \e, /). Using Lemma 7 we have 
that p is separable iff p is separable. But p is supported 
on (D 2 ® (F^ -1 and has rank N — 1, and therefore it is 
separable by assumption. □ 

Corollary 3a: If pis supported on <D 2 <&([! N and r(p) = 
N, then p can be written as a convex sum of projectors 
on N product vectors. 

Proof: Just given by the construction of the Proof of 
Theorem 1. 

Corollary 3b: If pis supported on (D 2 ®(D N and r(p) = 
N then r(p TA ) = N. 

Proof: This is a direct application of the previous corol- 
lary and Lemma 12. 

Corollary 3c: If p is supported on <D 2 ® <D N and is 
not separable, then r{p) > N. 

Proof: Direct application of Lemma 12 and Theorem 

1. 

B. rank(p)+rank(p T - 4 ) < 3iV 

We assume that r(p), r(p TA ) > N. We can use Lemma 
11 with Hi = R(p) and H 2 = R{p TA ) to determine the 
conditions under which there exist product vectors that 
can be subtracted from p. The problem reduces to finding 



the roots of a polynomial in a and a*. As explained in 
Section III, in the generic case these roots are also roots 
of another polynomial in a. If the corresponding projec- 
tors Pi are linearly independent [which imposes that the 
number of such product vectors be smaller than the min- 
imum of r{p) 2 ) and r( ( o T - 4 ) 2 ], then Lemma 11 provides 
us with a necessary and sufficient condition for separa- 
bility. We will illustrate this with the case <Z7 2 <g)(F 4 in the 
following section. 

Let us now show that if N < 10 then for a generic p 
we can always find less than the minimum of r(p) 2 and 
r(p TA ) 2 product vectors \e, f) G R(p) and with \e*,f) G 
R{p). Without loss of generality we can take k(p) = 
X > k(p TA ) = Y (otherwise we can substitute a — > a* in 
what follows). We also assume that there are no product 
vectors in K{p) since otherwise we can trivially reduce 
the problem to (D 2 (gitf^ -1 decreasing the ranks of p and 
P Ta via Lemma 7. We therefore have N > Y > X > 
and X + Y > N. Let us analyze separately the cases: 
(i) X + Y = N and; (ii) X + Y > N. If X + Y = N 
we have < X < [N/2], where [W] denotes the integer 
part of W. According of Lemma 11, the problem reduces 
to finding roots of a polynomial of degree X in a and Y 
in a*. In Appendix A we give a procedure to find those 
roots, which states that the number of roots (for a generic 
case) will not exceed 2 X ~ 1 [N 2 + (N - X) 2 — X(N — X)]. 
This number should be smaller or equal to (N + X) 2 . 
One can check that for N < 10 this condition is always 
fulfilled, (ii) If X + Y > N then we have [(N + l)/2] < 

Y < N. According to Lemma 11 there will be two (or 
more) polynomial equations of degree N — Y in a and 

Y in a*. In Appendix B we show that in that case the 
number of roots cannot exceed 2 N ~ Y Y. This number 
should be smaller or equal to (2N — Y) 2 , which is always 
the case for N < 10. 

Finally, using the results of the previous subsection we 
can ensure that if r(p) = r(p TA ) = N + 1 then either 
the polynomials of Lemma 11 have exactly N + 1 so- 
lutions, in which case p is separable and expresable in 
terms of TV + 1 product vectors, or it has none and then 
p is bound entangled. This is so because if one can sub- 
tract one product vector, then one recovers the situation 
of Corollary 3a. 

C. rank( ( o)+rank(p TA ) > 3iV 

We can use Lemma 11 with Hi = R(p) and H 2 — 
R(p TA ) to show that there always exist product vectors 
that can be subtracted from p. We can do that until we 
reach r(p) + r{p TA ) < 3 A. In that case, we can use the 
criterion given in the previous subsection to find out if the 
remaining operator is separable. Thus, this procedure 
gives us a constructive sufficient criterion for separability. 
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D. Invariance under partial transposition: p — p T 

Theorem 2: If p = p TA then p is separable. 

Proof: We prove it by induction. First, if p is sup- 
ported on <E 2 Cg) (D 1 then it is obviously true. Now, let us 
assume that it is valid if p is supported on(D 2 <S>^J N ~ 1 and 
let us prove that then it is also valid if it is on (D 2 ®(C N . 
According to Lemma 12, we just have to consider two 
cases: (i) r(p) — N: Using Theorem 1 we have that p is 
separable, (ii) r{p) > N: according to Corollary 2 there 
is a product vector \e T ,g) with |e r ) = |e*). Using Corol- 
lary 1 we can use this product vector to reduce the rank 
of p. But since \e T , g) = \e*,g) we have that the resulting 
operator is also equal to its partial transpose. We can 
proceed in this way until r(p) = N, in which case we use 
Theorem 1 to show that p is separable. □ 

Note that we could have chosen a different basis for 
the partial transposition. Moreover, since the property 
of separability is independent of invcrtible local transfor- 
mations A acting on (D 2 , we have the following 

Corollary 4: If [(A®I) p(A®I)i] TA = {A®I)p(A®I)\ 
for some non singular operator A, then p is separable. 

Theorem 2 suggests that if p is not very different 
from p T2 , then it should also be separable. Indeed, one 
can construct a powerful sufficient separability condition 
based on that theorem. For that we have to introduce 
some definitions. In 2 x N we can always write p as 

P=^ + ^ = p s + <^, (22) 

Here 2p s = P + P Ta , = i(\0) A (l\ - |1) A <0|), and 
AB = 4B t = tr A [o-y(p - p TA )]. This operator B can be 
decomposed as 

K 

B = J2*ih)(vil (23) 

i=l 

In particular, one of such is the spectral decomposition. 
Given one of such decomposition {A^, \vi)}fL 1 and a set 
of real numbers {ai}fL 1 we define the operator 

K 

C(a,X,v) = ]T |Ai|(a 2 |0)(0| + ar 2 \l)(l\) ® \v t )(v t \, (24) 

2 = 1 

which is obviously positive. We have: 

Theorem 3: Given a decomposition of B 
{Aj, \vi)}f =1 and a set of real numbers {ai}f =11 if 
\\C l l 2 {a, A, v)p^ 1 ^ 2 \\ 2 < 1, then p is separable. 

Proof: We define p s = p s — C(a, A, v) = p T s A > ac- 
cording to Lemma 1. Using Theorem 2 we have that p s 
is separable. Let |iWj) = at\0) — ia^~ sign(Aj)|l). Then, it 
is easy to check that 

K 

P = Ps + /,\^i\\wi,Vi)(wi,Vi\. (25) 

i=l 



which shows that p is separable. □ 

Thus, we can show that a density operator is separable 
if we can find a decomposition of B and a set of real 
numbers that fulfill certain conditions. In particular we 
can take the spectral decomposition of B and <ij = 1. 
Using the fact that \\AB\\ < \\A\\\\B\\ one can easily prove 
the following: 

Corollary 5: If p + p TA is of full range and \\(p + 
P Ta )~ 1 II II/ 0- P Ta II -• 1, then p is separable. This corollary 
implies that if p is full range and is very close to p TA then 
it is separable. 

Finally, let us mention that we have taken here all par- 
tial transposes with respect to the first system. One may 
wonder whether if p is invariant under partial transposi- 
tion with respect to the second system, then Theorem 2 
still holds. Actually, this is not the case. A counterex- 
ample can be found in Ref . (§) , by simply noting that the 
bound entangled state given there is equal to its partial 
transpose if we take that in a different basis. 

V. EXAMPLE: (F 2 ®<T 4 

As an illustration of the methods developed in the pre- 
vious sections to deal with separability when one has low 
rank density operators, here we will show how one can 
apply them in the case of (D 2 ®(C A . We will assume that: 
(i) p has no product vector in its kernel since otherwise 
using Lemma 7 we can reduce the problem to (U 2 (8>(F 3 and 
there we know that p is separable; (ii) r(p),r(p TA ) > 4 
since otherwise using Theorem 1 we know that p is sep- 
arable; (iii) r(p) + r(p TA ) < 12 since in that case we can 
use our separability criterion. This section is divided into 
3 parts: 

• r(p) — r(p TA ) = 5 

• r(p) + r(p TA ) < 12 and r(p) ^ r(p TA ). 

• r(p) — r(p TA ) = 6 

As before, we will denote by i = 1 . . .k(p)} and 

{|* 2 ),i = l...k(p TA )} a basis in K{p) and K(p TA ), 
respectively. We will use the fact that, according to 
Lemma 11, there exists a product vector |e, /) £ R(p) 
and \e*,f) £ R(p TA ) iff we can find an a such that 
there are at most 3 linearly independent vectors among 
{a(*lJO) + K |l>,«*(v&2jo) + 

A. rank(p)= =rank(p TA )= 5 

Each of the two kernels has dimension 3 and therefore, 
according to Lemma 11, in order to find subtractable 
product vectors we have to solve M(a, a*)f = 0, where 
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M(a,a*) 



( a(*i|0> - 
a(*a|0)- 

«<*3|0) " 

«*(*i|0> 
a*(*5|0> 
V a*<*||0) 



(*J|i> \ 
<**|i> 

-<*?|i> 



Hi) 



(*§ii>y 



(26) 



is a 6 x 4 matrix. Using Lemma 9 we know that the last 
three rows in M are linearly independent. Thus, this 
equation has non trivial solutions iff 

deti(a, a*) = det[Mi(a, a*)} = (27a) 
det 2 (a, a*) = det[M 2 (a, a*)] = (27b) 
det 3 (a, a*) = det[M 3 (a, a*)] = 0, (27c) 

where Mi. 2 , 3 are 4x4 matrices formed by the last three 
rows of M and the 1-st, 2-nd, and 3-rd rows, respec- 
tively. Note that deti, det 2 and det3 are polynomials 
linear in a and cubic in a* . Thus, we can write them as 



deti = aP£(a*) + P 3 ° (a*) = 
det 2 = aP£(a*) + P 3 °(a*) = 
det 3 = aP£(a*) + P 3 («*) = 0, 



(28a) 
(28b) 
(28c) 



where the subscripts denote the degree of the polynomi- 
als. In order to solve this equation we treat a and a* 



as two independent variables. Multiplying Eq. (28a) by 
P 3 (a*), Eq. (28b) by P 3 1 (a*) and subtracting them we 
find an polynomial equation for a* of degree 6. Doing the 



same thing but with Eq. (28c) instead of Eq. (28b) gives 
another polynomial equation for a* of 6th degree. Sub- 
tracting these two polynomials we obtain a polynomial 
of degree 5. Thus, there are at most 5 solutions (i.e. 5 
product vectors, as stated in Section IV B). Note that the 
solutions have still to fulfill all Eq. (^). Note also that 
if we know a we can calculate the kernel of the matrix 
A for this particular value and find so the corresponding 

If)- 



B. r(p)+r(p TA )< 12 where r(p) / r(p TA ) 

We are going to divide this part itself into: 

• r(p)= 6 and r(p TA )— 5 or r(p)= 5 and r(p TA )= 6, 

• r(p)= 7 and r(p TA )= 5 or r(p)= 5 and i(p TA )= 7. 

Let us begin with the case where r(p)= 6 and r(p TA )= 5. 
Then the dimension of the kernel of p is 2 and the one of 
the kernel of p TA is 3. Thus we can write the equations 
in the same way as in the previous case, but there are 
only two of them. We have 



deti = det[Mi(a,a*)] = 
det 2 = det[M 2 (a,a*)] = 0, 



(29a) 
(29b) 



where Mi and M 2 are the same matrices as above. Thus 
the way of calculating the a's is exactly the same as be- 
fore, obtaining at most 6 solutions for a*. On the other 
hand it is clear that the same can be done in the case 
where r(/?)=5 and r(p TA )— 6. 

Now we consider the case r(p)= 7 and i(p TA )— 5. Here 
we have only one equation that has to be fulfilled. We 
find 



deti = det[Mi(a,a*)] = 0, 
which leads to 



= aP} 



P 3 V 



(30) 



(31a) 



= {a*fQ\(a) + {a*) 2 Ql{a) + a*Q\(a) + Q?(a) (31b) 



Taking the complex conjugate of Eq. ( |3 1 a| ) we have: 

«*P 3 1 (a) + P 3 °(a) = 0, (32) 

where P is the same polynomial as P but with the coeffi- 
cients complex conjugated. As explained in the Appendix 
we are using now Eq. ([52]) to transform the polynomial 
in Eq. ( [31a| ) in one depending only on a. That is, we 
multiply Eq. ( pb| ) by P£(a) and Eq. @ by Q\{a) and 
subtract them. We do the same two more times and end 
up with a polynomial of 10-th degree in a. For the dif- 
ferent values of a we can now calculate the kernel of the 
matrix M and obtain |/}. Thus we have all the product 
vectors \e, f) fulfilling \e, f) G R(p) and \e*J) 6 R{p TA ). 

C. rank(p)=rank(p TA )= 6 

Here we have two vectors in both, the kernel of p and 
the one of p TA . Thus the equation that has to be fulfilled 
is 

= det(M(a, a*)) = a 2 P 2 2 (a*) + aP${a*) + P 2 °(a*) (33a) 
= {a*) 2 Ql{a) + a*Q\{a) + Q 2 {a). (33b) 



Again, we take the complex conjugate of Eq. (33a) and 
deal with the following two equations: 



(a*) 2 QlH + a*Ql(a) 
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(a*) 2 Ql{a) + a*Ql(a) + Ql{a) 



?§(«) 



(34a) 
(34b) 



First, we multiply Eq. ( paj ) by Q%(a), Eq. ( pb ) by 
Q 2 {a) and subtract them. Then we multiply Eq. (34a) 
by Q 2 ( a )i Eq- (34b) by Q\{a) and subtract them and 
divide the resulting polynomial by a*. Thus we find 



a*[Q 1 2 (a)Q 2 2 (a)-Ql(a)Q 1 2 (a)} + 

Q° 2 (a)Q 2 2 (a)-Q° 2 (a)Ql(a) = (35a) 
a*[Ql{a)Q%{a)-Ql{a)Q° 2 {a)] + 

Q 1 2 {a)Q a 2 { a )-Q Q 2 {a)Q 2 {a) = Q. (35b) 
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Now we multiply Eq. ( |35a[ ) by 
Q 2 2 (a)Q° 2 (a)}, Eq. (BEfl) by [Ql(a)Q|(a) 



[Ql(a)Q°(a) 



and subtract them again and end up with a polynomial 
of degree 8 in a. Note that this means that it may not be 
possible to write the density operator as a convex sum of 
6 product vectors. In fact, we have confirmed that and 
found examples in which one needs 8 product vectors. 



APPENDIX: REDUCTION OF ONE 
POLYNOMIAL 

In this appendix we show that given a generic poly- 
nomial Px,y(cx, a*) of degree X and Y in a and a* re- 
spectively, there are at most 2 X ~ 1 [A + Y(Y - X + 1)] 
roots, where without loss of generality we have assumed 
Y > X. The idea is to find another polynomial Q(a) 
of degree 2 X ~ 1 [X + Y(Y - X + 1)] such that all the 
roots of P are also roots of Q. Since Q has at most 
2 X ~ 1 [X + Y(Y - X + 1)] roots, this will prove our state- 
ment. 

We start by writing the equation P = as 

X 



= Px, Y (a,a*) = J2 akp Y( a *) ( Ala ) 

k=0 
Y 

= 5>*m(«), (Alb) 



fc=0 



where P Y and Q x denote polynomials of degree Y and 
X, respectively. We complex conjugate the last line in 
this equation and obtain 



= 



Y 

E< 

fc=0 



(A2) 



where Q is the same polynomial as Q but with the co- 
efficients complex conjugated. From now on we treat a 
and (3 = a* as t wo in depe ndent variables. Thus, we can 
write Eqs. ( Ala ) and (A2) as 



x 



J2^ k P Y ((3) = 

k=0 
Y 

J2^ k Q k xW = o. 



(A3a) 
(A3b) 



fe=0 



We proceed now in two steps: 

(1) If Y ^ X we use Eq. ( |A3a| ) to transform Eq. ( |A3b| ) 
to one of the form 



x 



]Ta fe i?!(/?) = 0, 



(A4) 



where Z = X+Y(Y — X). This can be done in the follow- 



ing way. First we multiply Eq. (A3a) by a' x Q x {f3), 



Eq. ( |A3b ) by Py(f3), and subtract them. In this way we 
obtain a polynomial equation of degree Y — 1 in a and 
Y + X in 0. We use again Eq. (A3a) to transform this 
one to a polynomial equation of degree Y — 2 in a. We 
proceed in the same way until we reach Eq. (A4). Note 
that if Y = X then we have automatically Eq. ( A4) with 
Z = X. 

(2) First, we multiply Eq. jA3a|) by R x (/3), Eq. (lAj) 



by Py {(3) and subtract them. In parallel, we multiply 



Eq. flA3a| ) by R° z ((3), Eq. (jA4j) by P Y {(3), subtract them, 
and divide by a the result. In this way, we obtain two 
polynomial equations of degree X — 1 in a and Y+Z in (3. 
We can do the same thing with these two, and proceed in 
the same vein until we reach a polynomial equation of the 
form Q(/3) = (i.e. Q(a) = 0) of degree 2^ X - 1 \Z + Y). 
All the roots of the original polynomial must be roots of 
Q. 



APPENDIX B: REDUCTION OF TWO 
POLYNOMIALS 

Following a similar procedure as in the previous ap- 
pendix, we show here that given two generic polynomials 
Px,Y and Px,y) of degree X and Y in a and a* respec- 
tively, there are at most 2 X Y common roots, where with- 
out loss of generality we have assumed Y > X. Treating 
a and j3 = a* as independent variables, if we write 



A" 



= P X , Y (a, a*) =J2 a k P Y (P), 



(Bla) 



X 



Q = P xx {a,a*)=Y,a k P Y (f3) 1 (Bib) 

(Blc) 



k=t) 



we can use the same procedure presented in step (2) in 
the previous appendix to obtain a polynomial equation 
Q(a) of degree 2 X Y. 

We emphasize that if one has more polynomials, one 
can further decrease the degree of Q(a). An example is 
shown in Section V, when the ranks of p and p TA equal 
5. 

This work has been supported by Deutsche Forshungs- 
gemeinschaft (SFB 407) and Schwcrpunkt "Quantcnin- 
formationsverarbeitung" , the Osterreichisher Fonds zur 
Frderung der wissenschaftlichen Forschung (SFB Pll), 
the European TMR network ERB-FMRX-CT96-0087, 
and the Institute for Quantum Information Gmbh. J. 
I. C. thanks the University of Hannover for hospitality. 
We thank, P. Horodecki, A. Sanpera and G. Vidal for 
discussions. 



10 



[1] R. F. Werner Phys. Rev. A 40, 4277 (1989). 

[2] See, for example, A. Peres, Quantum Theory: concepts 

and methods, (Kluwer Academic Publisheres, Dordrecht, 

1993). 

[3] See, for example, D. P. DiVincenzo, Science 270,255 
(1995). 

[4] R. Horodecki, P. Horodecki, and M. Horodecki, Phys. 

Lett. A230, 377 (1996). 
[5] A. Peres, Phys. Rev. Lett 77, 1413 (1996). 
[6] M. Horodecki, P. Horodecki and R. Horodecki, Phys. 

Lett. A223, 1 (1996). 
[7] E. Stromer, Acta Math. 110, 233 (1963); S. L. Woronow- 

icz, Rep. Math. Phys. 10, 165 (1976). 
[8] P. Horodecki, Phys.Lett. A232, 333 (1997); M. 

Horodecki, P. Horodecki, R. Horodecki, Phys. Rev. Lett. 

80, 5239 (1998). 
[9] C. H. Bennett, D. P. DiVincenzo, T. Mor, P. W. Shor, 

J. A. Smolin, and B. M. Terhal, Phys. Rev. Lett. 83, 

3081 (1999); D. P. DiVin cenzo, T. Mor, P. W . Shor, J. 



A. Smolin, B. M. Terhal, 



quant- ph/9908070 B. Terhal, 



quant-ph/9810091; C. H. Bennett, D. P. DiVincenzo, Ch. 



A. Fuchs, T. Mo r, E. Rains, P. W. Shor, J. A. Smolin, 



W. K. Wootters, fcuant-ph/9804053[ 
[10] R. Horodeck i, M. Horodecki, and P. Horodecki, [quant- 



ph/9811004. 



[11] M. Horodecki, and P. Horodecki, Phys. Rev. A59, 4206 
(1999). 

[12] K. Zyczkowski, P. Horodecki, A. Sanpera, and M. Lewen- 
stein, Phys. Rev. A58, 883 (1998); G. Vidal and R. Tar- 
rach, Phys. Rev. A59, 141 (1999); S. L. Braunstein, C. 
M. Caves, R. Jozsa, N. Linden, S. Popescu, R. Schack, 
Phys. Rev. Lett. 83, 1054 (1999). 

[13] P. H orodecki, J. A. Smo lin, B. M. Terhal, A. V. Thap- 
liyal, |quant-ph/9910122 . 

[14] W. Diir, J. I. Cirac, and R. Tarrach, Phys. Rev. Lett. 83 
3562 (1999). 

[15] M. Lewenstein and A. Sanpera, Phys. Rev. Lett. 80, 2261 
(1998); A. Sanpera, R. Tarrach, and G. Vidal, Phys. Rev. 
A58, 826 (1998), G. Vidal, PhD Thesis, Universitat de 
Barcelona (1999). 



11 



